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Abstract 

We present a new and simple method of calculating the occupation prob- 
ability of the number of total harmonic-oscillator quanta for a microscopic 
cluster-model wave function. Examples of applications are given to the recent 
calculations including a -|- n -|- n-model for ^He, a-|-t-|-n-|-n-model for ^Li, and 
a -|- a -|- n-model for ^Be as well as the classical calculations of a n-model 
for ^Li and a -\- a -\- a-model for ^^C. The analysis is found to be useful for 
quantifying the amount of excitations across the major shell as well as the 
degree of clustering. The origin of the antistretching effect is discussed. 

PACS number (s): 21.60.Gx, 23.40.Hc 
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Typeset using REVT^ 



The microscopic cluster model (MCM) is a many-nucleon theory which provides a unified 
picture of bound-state properties of nuclei and nuclear reactions. (See, for example, It is 
based on the assumption that the nucleons in the nuclei form substructures, called clusters. 
Though the MCM is capable of describing a variety of structure, its application has mostly 
been limited to two- or three-cluster system. Recent advances in the MCM have, however, 
enabled one to treat systems containing more than three clusters and thereby give a detailed 
description of light nuclei including halo nuclei . This extension of the applicability has 
been made possible by the inclusion of clusters other than the a-particle and by the use of 
the stochastic variational method The MCM, as a microscopic theory, attempts to 

derive the nuclear properties from a many-body Hamiltonian using a fully antisymmetrized 
wave function. It has, therefore, some relationship to the shell model. Efforts have been 
made to relate the wave function of the MCM to that of the SU{3) [0 or symplectic 
shell model. Such efforts are useful because one can compare the model space of the shell 
model and the MCM. It is hard to analyse a general MCM wave function in terms of shell- 
model configurations. We will show instead that it is easy to calculate the percentage of the 
harmonic-oscillator (HO) excitations involved in the MCM wave function. 

The MCM wave function consists of the antisymmetrized product of the intrinsic wave 
functions of the clusters and the relative motion functions between the clusters. The intrinsic 
wave functions are usually assumed to be described with simple configurations, whereas 
the relative motion functions are treated fiexibly enough to accommodate various types of 
correlations between the clusters as well as large spatial extension if necessary. The MCM 
wave function thus may contain a large number of HO quanta. To link it with the HO 
shell-model basis, it is useful to calculate the occupation probability of the number of total 
HO quanta. The occupation probability Pq of a definite number of total HO quanta Q for 
74-nucleon system is obtained by calculating the expectation value of the operator O 

= 1- f\e exp(^^ {Y.P^ [H^oii) - |] - Q)). (1) 

ZiTX Jo i=1 

Here Pi projects out either proton or neutron. It is set the unit operator when one calculates 
the number of total quanta occupied by both protons and neutrons. H}io{i) is the 3- 
dimensional HO Hamiltonian divided by hu = —7. 

The MCM wave function is conveniently generated from the Slater determinants of the 
Gaussian wave-packet single-particle (sp) functions 



with 




Here A is the antisymmetrizer and n = (ctiTi, a^T^) stands for the set of the spin- 
isospin quantum numbers of the nucleons. The Sj parameter or "generator" coordinate is 
a variational parameter in the generator coordinate method calculations or it is used in 
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an integral transformation [^,0 to derive the matrix elements between the Gaussian basis 



functions [10 



The matrix element of the operator O between the Slater determinants is given by 

(04si,...,s^)|O|0,,(s'i,...,s:,)) = ^ r dO eM-^QO) det{5}, (4) 

/7r Jo 

where the matrix B is defined by [i, j=l,...,y4) 

B,, = {^l,Jexv(jOP[Hno - (5) 

Since the constants, v and 7, are in general different, the calculation of the sp matrix element 
of Eq. (5) may seem difficult at first sight but can be performed rather easily with the use 
of the following formulae 



exp(ie[H„o - |)W(r) = exp(-2(l - z^)^) ^J.(r), (7) 

where z = e^^ . One can prove Eq. (7) by noting that the function, ^p1{v), is the generating 
function for 3-dimensional HO functions. Using Eqs. (6) and (7) in Eq. (5) yields the 
needed matrix element 



/ + 7 + (z/ - 7)2;2 2x , 4z/V A 

^ ^"Pp^ (^ + 7)^-(^-7)^^^ ^' ^ ^ (. + 7P-(^-7)^.-^ ' ■ ^ j^-'^-'' 

where z = 2; or 1 in accordance with {t\P\t) = 1 or 0. 

The value of v is usually chosen to give an appropriate size for the cluster, while the value 
of 7 is determined by the size of the whole nucleus. Hence the value of v is usually larger 
than that of 7. The integral in Eq. (6) then does not converge. However, even in this case 
Eq. (8) may safely be used because the matrix element should be continuous at all values 
of V and 7 and should have a finite value. We have confirmed the validity of Eq. (8) by 
comparing the occupation probability calculated by two ways for the Gaussian wave-packet, 
yjg: One is to use Eq. (8) and to change the integration variable 6 to the complex variable 
2; on a unit circle \z\ = 1. The other is to calculate directly the overlap between ip'^ and 
3-dimensional HO functions. Both methods have given the same result. 

The sp matrix element of i/no itself 

xe-t(-«')^ 5,,,, {r\P\r') (9) 
is needed to calculate the average number of total HO quanta contained in the wave function. 



Recently this quantity is used in Ref. ITl 
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The summation in the exponent of Eq. (1) runs over all the nucleons and the probability 
calculated with it in general contains the contribution from the center-of-mass (cm) motion 
unless the wave function is free from the spurious cm motion. Though the Slater deter- 
minants of Eq. (2) contain the dependence on the cm variable, our MCM wave functions 
generated from them by an integral transformation turn out to be free from any problem 
with the spurious cm motion [||. The probability calculated below is thus a purely intrinsic 
quantity. 

A generalization to a combined occupation probability is straightforward. For example, 
the probability -Pqi,q2 ^^at protons have Qi quanta and neutrons Q2 quanta or spin-up 
nucleons have Qi quanta and spin-down nucleons Q2 quanta is obtained by using Eq. (1) 
twice and noting the commutability of the corresponding operators. 



As an illustrative example let us consider Brink's a + a-model for Be [|T2[. The intrinsic 
wave function of the a-particle is constructed from the Slater determinant of a Os HO 
function with a size parameter u. When the two a-particles are separated by S and their 
relative orbital angular momentum is L, Pq is calculated by 



with 



4p N 12 



( l+p+(l- p)Z^ X 

/ • / %pdz Apdz , „ ^ 



where p = ■y /u, d = uS'^, and iL{z) = y (z) are the Modified Spherical Bessel functions 

of the first kind. Since f{z) is analytic in the unit circle, Cauchy's integral formula can be 
applied in Eq. (10) to yield Pq = f^'^\0)/Q\. The function f{z) has a leading term 
proportional to z^^^^^'^^ near z = 0. Hence Pq vanishes for Q < max(L,4), which is the 
consequence of the Pauli principle. Figure 1 shows the Pq values (Q=4, 6, and 8) for L=0, 
2, 4, and 6 as a function of S, with a choice of u = 0.25 fm~^ and 7 = 0.15 fm~^. The 
Q-dependence is not the same within the L=0, 2, and 4 members but different at each 
S particularly in the interval of 2~4 fm, where the energy functional is known to have a 
minimum. A maximum of Pq=4 appears around 5*^2.9 fm in case of L=0, while it shifts 
to a smaller separation of S^2.2 fm in case of L=4. The diagonal energy curve of ^Be as 
a function of 5* is expected to have a local minimum around the point where Pq=4 reaches 
a maximum. Then the difference in the Q-dependence on L and S values accounts for 



the antistretching |]13l that the energy functional of the cluster-model calculation shows a 
shrinkage of the cluster separation when L approaches the value of a band termination. 

Table I lists the Pq values in % for nucleons, protons, and neutrons for some of the wave 
functions obtained in our recent MCM calculations PJ5|JT^ using the Minnesota potential 
1511. A common value of z/=0.26 fm~^ is used to describe the intrinsic wave functions of 



a,t, and h clusters. The choice of 7 has some influence on the probability. It is set 0.17 
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fm~^ {huj = 14.4 MeV), a standard value used in a shell-model calculation for p-shell nuclei. 
The value is, however, reduced to 0.13 and 0.14 fm~^ for ^He and ^Li, respectively, to take 
into account their significantly large sizes. For the sake of reference the calculated root- 
mean-square (rms) radii for nucleons (matters), protons, and neutrons are included in the 
table. The Pq values are given as a function of Qexc = Q — Qmin, where Qmin is the minimum 
number of HO quanta for the lowest Pauli-allowed configuration. The lowest Ohuj component 
is around 50-60 % for most cases and the sum of 0, 2, and Ahuj components accumulates 
to about 90 %. The admixtures of higher components than Qexc = 4 are significant in the 
ground states of ^Li and ^Be and even larger in the ground state of ^He corresponding to 
its extended halo structure [4(b)]. The probability distribution spreads out to a very large 
number of HO quanta in ^Be and the 0^ state of ^^C, well-known cluster states. They are 
described as a bound state in a large basis. Our wave functions for ^^C are similar to those 
of Ref. |jl6|, which reproduces many properties of ^^C in the 3a-model. The parameter 



of the Minnesota potential is set to reproduce the energy of the O2 state. The ground 
state energy becomes then about 4.5 MeV lower than experiment. The calculated monopole 
matrix element is 4.0 fm^, which is in reasonable agreement with the experimental value of 
5.4±0.2 fm^. It is noted that no component is dominant in the 0^ state of ^^C. Of course 
it would be possible to maximize the probability with lower Q by choosing a smaller value 
of 7. However, the probability distribution would then spread to higher HO quanta in the 
ground state of ^^C. A recent large-basis shell-model calculation fl^ also suggests that the 
low-lying states in light nuclei from ^He to ''Li have significant admixtures of large HO 
quanta. Since the calculation seems promising, it is interesting to see how the 02" state of 
^^C is incorporated in its framework. It is also noted that the components with odd Qexc 
values for protons or neutrons are generally smaller. For example, in ^He -PQ<,xc(p)=i,Qcxc(n)=i 
is about 3 %, whereas Po,2 and P2,o are 9 and 11 %, respectively, and, among the probability 
of 16 % for Ahu excitations, the probability with Qcxdp) {Qexc{n)) =1 or 3 is less than 4 %. 

Comparing the results for ^~^Li, we see that the probability for neutrons has a larger 
change in the isotopes than that for protons. The change follows that of the neutron radius, 
which is consistent with the change of the neutron separation energy. In fact the nucleus ®Li 
has the smallest neutron separation energy among the three. Since the MCM consistently 



predicts the largest neutron rms radius for Li its average value, < Qexc >, for neutrons 



is largest among the three. A comparison of our result with that of Ref. |]TT| indicates that 
the latter wave functions, giving generally much smaller < Qexc > values, are rather close 
to simple shell-model configurations; e.g., < Qexc > for neutrons in their paper is about 0.1 
for ^Li and 0.5 even for ^Be. This may not be surprising because they use a single Slater 
determinant of the Gaussian wave-packets sp functions. 

The calculational method developed here has nothing to do with the assumption of the 
existence of clusters and can obviously be applied to those precise wave functions for few- 
nucleon system which are obtained with a sophisticated technique. For example, it will be 
interesting to analyse some of the wave functions obtained in Ref. 0. 

In summary, we have presented a new and simple method of calculating the occupation 
probability of the number of harmonic-oscillator quanta. It has been applied to the analysis 
of some of the wave functions obtained in a microscopic multicluster-model calculation. The 
analysis is found to be useful for clarifying the difference from the shell-model truncation 
made on the basis of a major shell as well as quantifying the degree of clustering or spa- 
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tial extension. The origin of the antistretching effect is discussed from the change of the 
occupation probabihty in the rotational band. 
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TABLES 

TABLE I. The occupation probability of the number of harmonic-oscillator quanta for micro- 
scopic multiclustcr-model wave functions. The probabilities for nucleons, protons, and neutrons 
are given in % in the upper, middle, and lower rows, respectively, as a function of oscillator exci- 
tations. When the probability for neutrons is the same as that for protons, only the proton case 
is shown. Asterisk indicates the probability of less than 1 % and dashed line represents vanishing 
probability. The average number of oscillator excitations is given in the column labeled < Qexc >• 
The details of the wave functions are referred to Ref. [4(b)] for ^He and ^Li, to Ref. [5] for ^Li, 
^Li, ^Li, and ^C, and to Ref. [14] for ^Be. 
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Figure Captions 



Fig. 1. The occupation probability of the number of oscillator quanta Q for the L = — 6 
states of ^Be. The wave functions of ^Be are assumed to be given by Brink's 2Q!-model 
[12] with a mean separation S. 
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